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Verfahren zum Visualisieren eines raumlich aufaelosten Datensatzes unter 
Verwenduna eines BeleuchtunQsmodells 

Die Erfindung betrifft ein Verfahren zum Visualisieren eines raumlich 
10 aufgeidsten Datensatzes unter Verwendung eines Beieuchtungsmbdells, 

sowie die VenA^endung dieses Verfahrens zum Generieren dreidimensionaler 
Darstellungen eines Korpers gemass dem Oberbegriff des unabhangigen 
Anspruchs der jeweiligen Kategorie. 

Zur Darstellung von drei- und hdherdimensionalen Datensatzen bedient man 
15 sich haufig zweidimensionaler Projektionen dieser Volumina oder 

Hypervolumina, welche dann graphisch ausgegeben und vom Betrachter 
interpretiert werden konneh. Fur den durchaus haufigen Fall, dass das 
Volumen in karthesicher bzw. isotroper Form vorliegt, und damit die Volumen 
in einem orthogonalen Raster angeordnet sind, das in alien drei 
20 Raumrichtungen dieselbe Auflosung besitzt, verfugen viele 

Grafikbeschleuniger modemer Datenverarbeitungsanlagen uber Hardware, 
die sich unmittelbar effizient nutzen iasst, um die Volumenvisuaiisierung 
durchzufuhren. Eine solche Ausstattung mit Grafikbeschleunigem in Form von 
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Grafikkarten ist auch bei handelsublichen Personalcomputern inzwischen 
durchgangiger Standard. 

Die Visualisierung raumlich aufgeloster Datensatze durch dreidimensionale 
bildliche Darstellungen, beispielsweise in einer Ebene, gewinnt dabei auf 

5 vielen technischen Gebleten zunehmend an Bedeutung. Das betrifft sowohl 
Animationen, beispielsweise bei Computerspielen Oder in der Werbung, als 
auch den industriellen Bereich und insbesondere die modeme medizinische 
Diagnostik und Therapie. Hier sind zahlreiche bildgeben 
Untersuchungsmetiioden, wie unter anderen die Computertomographie. die 

10 Kemspintomographie Oder Verfahren der Ultraschalltechnik bekannt, die, 
anhand von messtechnisch erfassten DatensStzen, Darstellungen von 
bestimmten Bereichen des menschlichen Korpers, von Organen; dam inneren 
von Blutgefassen Oder dem Herzen, des menschlichen Schadels usw. geben 
sollen. Dabei kommt es zunehmend darauf an, sowohl im wesentlichen 

15 statische Abbildungen, als auch bewegte Vorgange moglichst in Echtzeit 
abzubilden. Das ist zum Beispiel bei der Beobachtung der Bewegung des 
Herzens mittels eines Katheters oder wenn eine entsprechende Messsonde, 
Z.B. eine Ultraschallsonde, bei einer Operation das Auge des Arztes ersetzen 
muss, von zentraler Bedeutung. Verwandte bildgebende Verfahren sind 

20 naturlich auch aus der industriellen Technik, beispielsweise zur 

zerstorungsfreien Uberprufung sicherheitsrelevanter Komponenten wie Rad- 
und Achssatze be! Schienenfahrzeugen, Druckbehaltern, Rohren und dunnen 
Leitungen, z.B. in der Kraftwerkstechnik und in vielen anderen Bereichen 
mehr wohlbekannt. 

25 Daruber hinaus gewinnt die Visualisierung raumlich aufgeloster Datensatze 
durch dreidimensionale bildliche Darstellungen, insbesondere von 
Radardaten, Sonardaten (Ortung und Navigation), von seismischen 
Datensatzen, von Wetterdaten oder beispielsweise die Visualisierung im 
Rahmen von Finite Elemente Analysen eine immer grossere Bedeutung. Aber 



auch fur Computersimulationen auf den verschiedensten Gebieten, unter 
anderem im Bereich der Radar- Ultraschall- und Sonartechnik gibt es 
zahllose Anwendungen. 

Zunehmend wird auch auf eine immer realistischere Darstellung der 
5 messtechnisch erfassten und dann projizierten Datensatze Wert gelegt. Das 
heisst, dass die Tendenz hin zu immer hoherer raumlicher Auflosung des zu 
beobachtenden Objekts geht, wobei es darauf ankommt, den durch eine 
Messapparatur erfassten, in der Regel dreidimensionaien Datensatz 
perspektivisch korrekt, beispielsweise auf die zweidimensionale Oberfldche 

10 eines Computermonitors zu projizieren, und dem projizierten Volumen in 
realistischer Weise Lichtreflexe einzuschreiben, um den menschlichen Visus 
be! der Orientierung in der rSumlichen Tiefe in einem dreidimensionaien 
graphischen Datensatz zu unterstutzen. Insbesondere bei stereoskopisclien 
Projektionen. wenn also dem mensclilichen Visus durch Uberiagerung von 

15 zwei in der Perspektive leicht gedrehte Projektionen durch eine geeignete 

Projektionseinrichtung ein dreidimensionales Bild vermittelt werden soil, ist die 
Einbeschreibung von Beleuchtungseffekten in die Projektion von besonderer 
Bedeutung. 

Dazu wird ein bestimmtes Beleuchtungsmodell zugrunde gelegt, dass durch 
20 sogenannte Beleuchit:ngsfunktionen beschrieben wird, die die 

Wechselwirkung von (virtuell) eingestrahltem Licht mit den Objekten im 
Volumen, einschliesslich DSmpfung, Reflexion und Streuung approximieren. 
Es liegt auf der Hand, dass dazu enorme Rechenleistungen erforderlich sind, 
die ohne weiteres auch nicht von den heutzutage auf dem Markt verfugbaren 
25 Rechnersystemen erbracht werden. 

Das aus dem Stand der Technik bekannte Prinzip der Visualisierung 
mehrdimensionaler graphischer Datensatze mit Hilfe handelstiblicher 
Graphikhardware unter Berucksichtigung eines bestimmten 



Beleuchtungsmodells soil im folgenden unter Bezugnahme auf die Fig. 1 bis 3 
kurz skizziert werden. Zur Unterscheidung des Stands der Technik von dem 
erfmdungsgemassen Verfahren sind in den Fig. 1 bis 3 die Bezugszeichen mit 
einem Hochkomma versehen. All den aus dem Stand der Technik bekannten 
Prinzipien zur Visuaiisierung von Datensatzen D'. die Daten D'(a', p', y') 
(gemessen in einem Messkoordinatensystem K'm mit Koordinatenachsen a'. 
P'. y' umfassen, ist gemeinsam. dass die Beleuchtungsfunktionen. die das 
Beleuchtungsmodell definieren. im Projektionsraum P' ausgewertet werden. 
da dort naturlich der Sehvektor S' des Beobachters bzw. der 
Beleuchtungsvektor definiert ist. Im Rahmen dieser Anmeldung wird der 
Vorgang des Berechnens der Beleuchtungsfunktionen in Aniehnung an die 
Nomenklatur der einschlagigen Literatur als "Shading" bezeichnet. Devon zu 
unterscheiden ist der im Rahmen dieser Anmeldung als "Rendern" 
bezeichnete Vorgang. Unter Rendern soli im folgenden die Linearisierung 
nach Schneiden In einem Originalvolumen bzw. in einem gemessenen 
Datensatz und das anschliessende Obersetzen der Schnittpunkte in die 
Geometrie eines Projektionsraumes verstanden werden. 

In einem einfachen Fall liegt, wie In Fig. 1 an einem bekannten Beispiel aus 
dem Stand der Technik schematisch dargestellt. der zu visualisierende 
Datensatz D', der durch eine entsprechende Vorschrift Oder eine 
Messvon-ichtung 1', wie belspielsweise durch einen Kemspintomographen 1' 
durch Messung in einem Originalvolumen Go' generiert wurde, in 
karthesischer oder isotroper Form vor. Das Originalvolumen Go' hat die Form 
eines Quaders und die Volumenelemente V sind in einem orthogonalen 
Raster angeordnet, das in alien drei Raumachsen a', p', y' des 
Messkoordinatensystem Kn, dieselbe Auflosung besitzt. Zur Projektion eines 
solchen Datensatzes D', belspielsweise auf einem Beobachtungsmonitor 2', 
wird der Datensatzes D' in eine Datenverarbeitungsanlage 3' geladen. Diese 
verfQgt uber handelsubliche Graphikbeschleiniger, die Graphikhardware 4' 
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umfassen, welche sich effizient nutzen lasst, um die Volumenvisualisiemng 
auszufuhren. 

Hierfur findet haufig das Konzept der Textur Verwendung, wobei die Textur 
durch einen Datensatz definiert ist, aus dem zweidimensionale polygonale 
5 Flachen herauskopiert, also herausgelesen werden kennen, was der oben 
erwahnten graphischen Schneideoperation entspricht. 

Dazu wird, wie in Fig. 2 schematisch gezeigt, bei der Abbildung eines digitalen 
Bildes 5' aus dem Orginalvolumen Go', d.h. aus einem digitalen Datensatz D', 
ein beliebiges Polygon auf ein anderes Polygon Ep im Projektionsraum P' 

10 abgebildet, wobei das Polygon Tp' des Originalvolumens Go' geometrisch 
nicht ahnlich zu dem Polygon Ep im Projektionsraum P' zu sein braucht. Wird 
wie im vorliegenden Beispiel ein Originalvolumunen Go' durch einen 
dreidimensionaler Datensatz D', der aus einer oder mehreren 2D~ oder 3D- 
Texturen Tp' aufgebaut ist, reprasentiert, so wird ein solcher 

15 dreidimensionaler Datensatz D' haufig auch als 3D-Textur bezeichnet. 

Es versteht sich, dass in einem besonders einfachen Fall die Textur Tp' 
identisch mit dem Polygon aus der Schneideoperation sein kann. 

Das Originalvolumen Go', also der Datensatz D', wird beispielsweise 
senkrecht zu einer Betrachtungsrichtung B* eines (virtuellen) Betrachters 6\ in 

20 eine bestimmte 2^hl von Texturen Tp' unterteilt, deren Ecken dann in die 

Betrachtungsposition gedreht, perspektivisch entzerrt und dargestellt werden. 
Diese Darstellung des Originalvolumens Go' als zweidimensionale Projektion 
wird durch wiederholtes "Schneiden" des Originalvolumens G©', d.h. des 
dreidimensionalen Datensatzes D' senkrecht zur Betrachtungsrichtung 6' 

25 erreicht. Das Schneiden besteht dabei darin, dass die Werte aus dem 
Datensatz D' interpoliert und in eine Projektionsebene, die in dem 
vorliegenden Fall mit dem Projektionsraum P' identisch ist, projiziert werden. 
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Die Menge der Schnitte, die hier zweidimensionale Bilder sind, mussen dann 
durch eine Integrationsvorschrift zu einem einzigen Bild 5' zusammengefugt 
werden. Die Integrationsvorschrift kann dabei eine simple punktweise 
Summation sein. oder durch eine bestimmte Vorschrift. die oft auch als 

5 Beleuchtungsfunktion bezeichnet wird, im Rahmen eines sogenannten 
Beleuchtungsmodells erfolgen. Die Beleuchtungsfunktionen, bzw. die 
Beleuchtungsmodelle berucksichtigen also durch den Vorgang des Shadings 
approximativ die Wechselwirkung von (real oder virtuell) eingestrahltem Licht. 
einschliesslich Dampfung. Reflexion und Streuung, mit den Objekten des 

10 Originalvolumens Go*, also dem dreidimensionalen Datensatz D'. 

Soil eine stereoskopische Projektion en-eicht werden, muss das zuvor 
geschilderte Verfahren zusatzlich noch in Bezug auf mindestens einen 
zweiten Sehvektor, dessen Richtung vom Sehvektor S' leicht verschieden ist. 
durchgefuhrt und einer geeigneten stereoskopischen Projektionseinrichtung 
15 zugefQhrt werden. D.h.. fur jedes Auge muss entsprechend seiner Position 
eine Ansicht berechnet werden. 

Die Literatur der Volumenvisualisierung kennt dabei zahlreiche Ansatze, wie 
einem zu visualisierenden Datensatz D' optische Eigenschaften 
zugeschrieben werden konnen. Die bekannteste Klasse solcher Verfahren 
20 betrachtet die Intensitat als eine Fomn der optischen Dichte. so dass 

Dichteschwankungen Lichtstreuung und Reflexion verursachen und die Dichte 
selbst lichtabsorbierend wird. In der Literatur werden Algorithmen. die auf 
dieser Hypothese bemhen, pauschal als Gradientenrenderer bezeichnet. 

Naturlichenweise wertet man die Beleuchtungsfunktionen im Projektionsraum 
25 P' aus. da dort der Sehvektor S' und der Beleuchtungsvektor definiert sind. 
Das heisst, das Originalvolumen Go" (also der Datensatz D') wird zunachst in 
Ebenen geschnitten, d.h. in Texturen Tp' zerlegt. und in den Projektionsraum 
P' abgebildet, ein Vorgang der auch als Texturremapping in der Literatur 
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bezeichnet wird. Anschliessend werden im Projektionsraum P' die 
Beleuchtungsfunktionen ausgewertet, das heisst, auf die projizierten (also 
gedrehten und entzerrten) Texturen Tp' im Projektionsraum P' angewendet. 

Dieses aus dem Stand der Technik bekannte Verfahren ist einerseits deshalb 
5 nachteilig, weil der Sehvektor S' bzw. der Beleuchtungsvektor weder auf den 
ebenen Texturen Tp' des Originalvolumens Go' noch auf den projizierten 
Ebenen Ep im Projektionsraum P' einen konstanten Betrag haben; vielmehr 
ist der Betrag des Sehvektors, wie im linken Bild von Fig. 3 schematisch 
dargestellt, auf Kugelschalen K* konstant. Somit andert sich der Sehvektor S' 
10 und der Beleuchtungsvektor im allgemeinen uber einen Schnitt, d.h. uber eine 
ebene Textur Tp* bzw. uber eine Ebene Ep im Projektionsraum P'. 

Daher bilden die Ebenen, die durch Sehvektoren S' mit konstantem Betrag 
gekennzeichnet sind, im Originalvolumen Go konzentrische Kugelschalen K'. 
Weil somit betragsmassige Anderungen des Sehvektors S' nicht linear von 

15 den Texturkoordinaten der Texturen Tp', bzw. der projizierten Ebenen Ep 
abhdngen, konnen diese auch nicht auf der Graphikhardware 4' unmittelbar 
kodiert werden. Daher muss bei den aus dem Stand der Technik bekannten 
Verfahren die projizierte Ebene Ep im Projektionsraum P\ wie im rechten Bild 
der Fig. 3 schematisch dargestellt, in kleinere Flachen zeriegt werden, wobei 

20 zwischen den Eckpunkten der Teiifldchen linear interpoliert wird. 

Neben verschiedenen Nachteilen, wie Erhohung der benotigten Rechenzeit, 
liegt ein besonders gravierender Nachteil darin, dass durch die zuvor 
beschriebene Interpolation im Projektionsraum P* zwangslaufig 
Linearisierungsfehler auftauchen. 

25 Wahrend die zuvor erwahnten Linearisierungsfehler bei der Abbildung eines 
im wesentlichen karthesischen und isotropen Originalvolumens Go', also eines 
karthesischen Datensatzes D', sich in manchen Fallen noch in einem 
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vertretbaren Rahmen bewegen, werden die Linearisierungsfehler bei nicht 
karthesischen Datensatzen D' unvertretbar hoch, so dass entweder l^eine 
brauchbare Abbildung im Projektionsraum P' entsteht oder der 
Rechenaufwand stelgt so dramatisch an. dass Echtzeitabbildungen, wie sie 

5 beispielsweise in der medlzinischen Technik haufig unerlassiicli sind. mit den 
aus dem Stand der Technik bekannten Projektionsverfahren und der derzeitig 
verfugbaren Hardware nicht mehr moglich sind. Es versteht sich, dass durch 
eine zusatzliche Auswertung der Beleuchtungsfunktionen im Projektionsraum 
P- der Rechenaufwand nochmals dramatisch steigt und bei stereoskopischer 

10 Projektion noch weiter erhOht wird. 

Dabei liegen gerade bei vielen fCir die Praxis wichtigen Anwendungen die 
Datensatze D'. die im Originalvolumen Go' gewonnen wurden, eben nicht in 
karthesischen, d.h. in rechtwinkligen Koordinaten vor. Das hat seine Ursache 
in erster Linie in der Art der Datenaquisiticn. Typische Beispiele sind die 
15 (Rontgen) Computertomography oder spezielle Ultraschalltechniken, bei 
welchen die zu projizierenden Datensatze D" beispielsweise in 
Zylinderkoordinaten vorliegen. Bei der Venvendung sehr moderner 
Abtastsysteme ist somit die relativ lange Zeit. die zur Erstellung einer 
dreidimensionalen Abbildung im Proiektionsraum P' mit den bekannten 
20 Verfahren aus dem Stand der Technik benotigt wird. weniger durch die 
Datenaquisiticn als solche, also die Erstellung des Datensatzes D' an sich 
verursacht. sondem vielmehr durch den Vorgang der Visualisiemng des 
Datensatzes D*. So sind bereits sehr schnell abtastende Ultraschallsonden 
bekannt. mit denen ein interessierendes Originalvolumen Go' in mehreren 
25 Ebenen gleichzeitig und kreisffirmig abgetastet wird. Solche Ultraschallsonden 
umfassen beispielsweise mehrere, urn eine Achse schwenkbare oder 
rotierbare, linear oder allgemein in einem Array angeordnete 
Ultraschallwandler, von denen einige oder alle parallel betrieben werden. so 
dass ein gleichzeitiges zylindersymmetrisches Abtasten in mehreren Ebenen 
30 simultan ermSglicht wird. Sehr haufig sind solche schnell abtastenden 
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Untraschallsonden allein schon aus EffizienzgrCinden nicht in der Lage, das 
Volumen karthesisch zu vermessen. 

Die so gewonnen Datensatze D' liegen somit in einer Zylindersymmetrie vor 
und konnen nicht unmittelbar auf herkommlich verfugbarer Graphikhardware 

5 kodiert werden. Wenn eine schnelle Visualisieaing der Datensatze D' 
notwendig ist, urn beispielsweise eine Echtzeitdarstellung bei 
Videofirequenzen mit typischen Bildlaufzeiten von weniger als 1/25 Sekunde 
zu en-eichen, verbieten sich prinzipiell komplizierte und zeitaufwendige 
Rechenoperationen zur Aufbereitung der gemessenen Datensatze D' fur die 

10 Verarbeitung in der Graphikhardware 4'. Somit scheidet insbesondere eine 
aufwendige Koordinatentransformation in ein karthesisches 
Koordinatensystem aus. 

Zur Losung dieses Problems wird in der EP 1 059 612 A1 ein Verfahren zum 
Visualtsieren eines raumlich aufgelosten Datensatzes D' angegeben, dass es 

15 durch Vermeidung aufWendiger Teilschritte auf besonders effiziente Weise 
eriaubt, auch nicht karthesische Datensatze D* auf allgemein verfugbarer 
Graphikhardware 4\ wie sie beispielsweise in handelsublichen 
Personalcomputern zum Einsatz kommt, enorm schnell zu verarbeiten und: so 
dreidimensionale Darstellungen, selbst von bewegten Objekten, in Echtzeit, 

20 d.h. bei typischen Videofrequenzen, zu projizieren und darzustellen. Das 
vorgenannte Verfahren ist in der EP 1 059 612 A1, deren Inhalt hiermit in 
diese Anmeldung aufgenommen wird, eingehend beschrieben und braucht 
daher nicht mehr im Detail eriautert zu werden. 

Das in der EP 1 059 612 A1 vorgeschlagene Verfahren erlaubt es zwar auf 
25 besonders elegante und effiziente Weise mehrdimensionale Datensatze D\ 
die in beliebigen krummlinigen Koordinaten, z.B. in Zylinder- Oder 
Kugelsymmetrie vorliegen, ohne grossen Rechenaufwand und enorm schnell 
auf gewdhnlicher Graphikhardware unmittelbar zu kodieren. Das Problem der 
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schnellen und effizienten Auswertung der Beleuchtungsfunktionen bleibt 
jedoch unberuhrt. 

Ausgehend von diesem Stand der Technik ist es daher eine Aufgabe der 
Erfindung. ein Verfahren zum Visualisieren eines raumlich aufgelosten 
Datensatzes bereitzustellen, das es auf besonders effiziente Art eriaubt die 
Beleuchtungsfunktionen auszuwerten, wobei die benatigte Rechenzeit im 
Vergleich zu Verfahren aus dem Stand der Technik enorm reduziert wird. 

Die diese Aufgabe ISsenden Gegenstande der Erfindung sind durch die 
Merkmale des unabhSngigen Anspruchs der jeweiligen Kategorie 
gekennzeichnet. 

Die abhangigen Anspruche beziehen sich auf besonders vorteilhafte 
Ausfuhrungsfomnen der Erfindung. 

Erfindungsgemass wird somit ein Verfahren zum Visualisieren eines raumlich 
aufgelosten Datensatzes unter Verwendung eines Beleuchtungsmodells 
vorgeschlagen, wobei ein Datum des Datensatzes jeweils einem 
Volumenelement zugeordnet ist. dessen Lage durch Koordinaten in einem 
Messkoordinatensystem beschrieben wird. Die Daten werden als mindestens 
eine Textur in eine Grafikhardware geladen. urn in einem Projektionsraum 
eine bildliche Darstellung zu erzeugen. Dabei wird das Beleuchtungsmodell im 
Messkoordinatensystem ausgewertet. 

Wesentlich fur die Erfindung ist es somit, dass ein fur die Visualisierung des 
Datensatzes zugrunde gelegtes Beleuchtungsmodell bzw. die das 
Beleuchtungsmodell definierenden Beleuchtungsfunktionen im 
Messkoordinatensystem des Orginalvolumens ausgewertet werden. Das 
5 heisst. der Vorgangs des Shadings findet vollstandig im Originalvolumen statt 
und ist vom Vorgang des Renderns, wie er eingangs definiert wurde. 
vollstandig getrennt. 



- 1 1 - 

Bevorzugt werden dabei die Daten des Datensatzes, die in einer Messung, 
beispielsweise einer Messung mittels einer Ultraschallsonde am Herzen, 
generiert wurden, ohne Transformation aus dem Messkoordinatensystem in 
ein anderes Koordinatensystem, insbesondere ohne Transformation in ein 
karthesisches und / Oder isotropes Koordinatensystem, verarbeitet. Das 
heisst, be! Anwendung des erfindungsgemdssen Verfahrens, be! welchem die 
Beleuchtungsfunktionen im Messkoordinatensystem ausgewertet werden, 
werden die Daten durch elementare Rechenoperationen so aufbereitet, dass 
sie ohne Koordinatentransformation in die Grafikhardware gejaden und von 
dieser verarbeitet werden konnen. 

Bevorzugt kommt das erfindungsgemasse Verfahren in einer 
Ausfuhrungsvariante zum Einsatz, bei dem das Messkoordinatensystem ein 
nicht-karthesisches Messkoordinatensystem ist. Im speziellen kann das 
Messkoordinatensystem ein Zylinder- oder Kugelkoordinatensystem Oder ein 
15 anderes nicht-karthesisches Koordinatensystem sein. So konnen die Daten 
beispielsweise mittels einer sehr schnell abtastenden. rotierenden 
Ultraschallsonde generiert werden. Mit solchen Ultraschallsonden ist es 
moglich ein interessierendes Originalvolumen, beispielsweise das Innere 
eines menschilchen Herzens, in mehreren Ebenen gleichzeitig und 
20 kreisformig abzutasten. Solche Ultraschallsonden konnen mehrere, um eine 
Achse schwenkbare oder rotierbare, linear Oder allgemein in einem Array 
angeordnete Ultraschallwandler umfassen, von denen einige parallel 
betrieben werden, so dass ein gleichzeitiges zylindersymmetrisches Abtasten 
in mehreren Ebenen ermoglicht wird. Somit besitzt das 
25 Messkoordinatensystem, in dem die Daten des Datensatzes vorliegen, 
ebenfalls Zylindersymmetrie. 



5 



10 



Es versteht sich, dass das erfindungsgemasse Verfahren in einer speziellen 
Ausfuhrungsvariante auch vorteilhaft auf Daten eines karthesischen 
Datensatzes angewendet werden kann. Somit ist das erfindungsgemasse 
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Verfahren keineswegs auf nicht-karthesische Messkoordinatensysteme 
beschrankt. 

In einem fur die Praxis wiclitigen AusfCihrungsbeispiel des 
erfindungsgemassen Verfahrens, wird zwischen den Daten des Datensatzes 
5 im Messkoordinalensystem linear interpoliert. Dadurch kSnnen die Daten 
ohne DurchfQhrung einer Koordinatentransformation in ein kartiiesisches 
System in die Grafikhardware geladen und von dieser verarbeitet werden. 

Nicht zuletzt dadurch, dass die Daten des Messdatensatzes, die in einem 
krummlinigen, nicht-karhteslsclien Messkoordinatensystem vorliegen kdnnen, 

10 beispielweise in Zylinderkoordinaten Oder in Kugelkoordinaten, in einer fur die 
Praxis sehr wichtigien AusfOhrungsvariante des erfindungsgemassen 
Verfaiirens ohne vorherige Koordinatentransfomnation in die Grafikhardware 
geladen und von dieser verarbeitet werden kGnnen, konnen die Daten des 
Datensatzes auch nahe einer Singularitat ausgewertet werden. Seiche 

15 Singularitaten findet man beispielsweise in Zylinderkoordinaten, die durch 

eine Koordinate, die einem Radius entspricht, eine Winkelkoordinate und eine 
weitere Raumkoordinate beschrieben werden konnen, bei einem Wert des 
Radius von Null. Das heisst unter der Aussage bei "Radius Null liegt eine 
Singularitat vor" ist zu verstehen, dass die Daten in einem 

20 Messkoordinatensystem mit Zylindersymmetrie in der Nahe bzw. exakt bei 
Punkten mit einem Wert der Radius-Koordinate von Null nicht durch eine 
Koordinatentransformation auswertbar sind. Da jedoch in der oben 
beschriebenen AusfOhrungsvariante keine Koordinatentransformation 
stattfindet, kennen durch Anwendung dieser Ausfuhmngsvariante des 

25 erfindungsgemassen Verfahrens auch Daten des Datensatzes nahe einer 
Singularitat ausgewertet werden. Selbstverstandlich ist das 
erfindungsgemasse Verfahren grundsatzlich auch anwendbar, wenn die 
Daten des Datensatzes, insbesondere zur Weiterverarbeitung in der 
Grafikhardware, einer Koordinatentransformation unterzogen werden. 
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Insbesondere konnen die Daten des Datensatzes in einer weiteren Variante 
des erfindungsgemassen Verfahrens eine volumenaufgeloste Abtastung eines 
Korpers, zum Beispiel eines Tells ernes menschlichen Korpers reprasentieren, 
so dass die bildliche Darstellung eine dreidimensionale Volumendarstellung, 
5 insbesondere unter anderem auch eine semi-transparente Darstellung des 
Kdrpers ist Solche Darstellungen sind unter anderem dann von Vortell, wenn 
sich der Betrachter im Volumen des dargestellten Korpers orientieren muss 
und / Oder wenn der K6rper oder Teile des dargestellten Kdrpers gewissen 
Bewegungsablaufen unterliegt. So kann das Verfahren unter anderem zur 

10 Beobachtung eines schlagenden Herzens oder bei einer Operation am 
schlagenden Herzen. zum Beispiel unter Venvendung einer 
Ultraschallmesseinrichtung, besonders vorteilhaft eingesetzt werden. Da es 
durch Einsatz des erfindungsgemassen Verfahrens moglich ist, enorm schnell 
dreidimensionale Darstellungen mit Hilfe einer gewohnlichen Grafikhardware 

15 zu generieren, sind dreidimensionale Abbildungen von sehr hoher Auflosung 
in Echtzeit moglich, auch mit Videofrequenzen mit Bildlaufzeiten von 
typischerweise 1/25 Sekunden. 

Selbstverstandlich ist es durch Einsatz des erfindungsgemassen Verfahreris, 
bzw. durch Einsatz einer oder mehrerer, d.h. durch geeignete Kombination^ 

20 der zuvor geschilderten Ausfuhrungsvarianten auch moglich, die bildliche 
Darstellung als stereoskopische Projektion zu erzeugen. Insbesondere bei 
solchen stereoskopischen Projektionen, wenn also dem menschlichen Visus 
durch Oberlagerung von zwei in der Perspektive leicht verschobenen 
Projektionen durch eine geeignete Projektionseinrichtung ein 

25 dreidimensionales Bild vermittelt werden soil, ist die Einbeschreibung von 
Beleuchtungseffekten in die Projektion von besonderer Bedeutung. Da die 
verschiedenen Ausfuhrungsvarianten des erfindungsgemassen Verfahrens 
eine enorm schnelle Verarbeitung der Bilddaten, d.h. insbesondere eine 
enorm schnelle und effiziente Auswertung der Beleuchtungsfunktionen Im 

30 Messkoordinatensystem eriauben, sind selbst stereoskopische Projektionen, 
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eventuell unter Zuhilfenahme geeigneter Projektionsvorrichtungen. wie z.B. 
geeigneter 3D-Brillen, sogar von Bewegungsablaufen in Echtzeit und mit 
extrem hoher Auflosung moglich. 

Somit eignet sich das Verfahren in seinen verschiedenen 
Ausfuhrungsvarianten insbesondere zur Verwendung fur medizinische 
ZwGcke, zum schnellen Generieren dreidlmensionaier Darstellungen eines 
Korpers, insbesondere eines menschlichen Kerpers oder Teilen davon, 
anhand von messtechnisch gewpnnenen Daten. 

Es versteht sich, dass auch jede geeignete Kombination der zuvor 
dargestellten Ausfuhrungsvarianten vorteilhafl einsetzbar ist und darOber 
hinaus das Verfahren nicht nur fur medizinische Zwecke. sondem ganz 
allgemein auch in der industriellen Technik, z.B. zur Untersuchung schwer 
zuganglicher Bereiche einer Aniage, hervorragend geeignet ist. Insbesondere 
mussen die Daten des zu visualisierenden Datensatzes nicht unbedingt durch 
eine messtechnische Einrichtung, wie beispielsweise eine Ultraschallsonde 
generiert werden, sondem konnen beispielsweise auch durch eine 
mathematische Vorschrift, eine Simulation oder anders zur Verfugung gestellt 
werden. 

Bevor die Erfindung anhand der Zeichnung naher eriautert wlrd, sollen zum 
besseren Verstandnis des vorliegenden Anmeldetextes die fur die Erfindung 
wesentlichen mathematischen Grundlagen zur Visualsiemng eines raumlich 
aufgelosten, insbesondere eines nicht-karthesischen Datensatzes eriautert 
werden. 

Da bei einer Abbildung Ebenen des Originalvolumens, wegen der 
Perspektive, nicht in Ebenen sondem in Kugelschalen im Projektionsraum 
abgebildet werden, andem sich der Sehvektor und der Beleuchtungsvektor 
uber einen projizierten Schnitt. Weil diese Anderung nicht linear, bestenfalls 
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naherungseise von den Koordinaten der Textur abhangen, kann sie nicht als 
Texturoperation auf der Hardware codiert werden. Daher wird die Textur im in 
kleinere Teilflachen zerlegt und zwischen den Eckpunkten der Teilflachen wird 
linear Interpoliert. Die Eckpunkte dieser Teilflachen werden auch als Vertices 
5 bezeichnet. Bine solche Zerlegung in kleinere Teilflachen wird im Rahmen 
dieser Anmeldung als Tesselierung bezeichnet Dabei findet die Tesselierung 
im Originalvolumen statt. Das Originalvolumen wird in konzentrische 
Kugelschalen geschnitten, so dass der Betrag des Betrachtungsvektor auf 
einer gegeben Textur, die einer herausgeschnittenen Kugelschale entspricht, 
10 konstant bleibt, so dass der Betrachtungsvektor nicht entzerrt werden muss. 
Dadurch werden die ubitcherweise auftretenden Auflosungsverluste durch die 
Entzerroperationen minimiert. Die Geometrie wird also im ursprunglichen 
Originalvolumen linearisiert, in dem erfindungsgemass auch die 
Beleuchtungsfunktionen ausgewertet werden. 

15 

Der maximale Fehler, welcher durch eine Tesslierung mit Winkelinkrementen 
von kleinen Winkein a bei einem Betrachtungsradius r entsteht ist 

err ^Al- col ^^^f^r— (eq. 1). In Masseinheiten von Pixein, ergibt sich so 



20 fur r = 200 also etwa 6 Grad. Dieser Fehler ist in karthesischen Volumina 
teilweise tolerierbar - in krummlinigen Koordinaten, wie beispielsweise in 
Systeme mit spharischer oder zylindrischer Symmetrie, vollig unvertretbar 
gross. 

25 Einen Ausweg bildet das sogenannte „lnnere Shading" gemass der 

vorliegenden Erfindung, dessen mathematische Grundlagen im folgenden 
naher eriSutert werden. Beim inneren Shading wird die Geometrie des 
Originalvolumens in ihren nicht-kartesischen Koordinaten linearisiert und die 
Vertices in die Geometrie des Projektionsraumes ubersetzt. Die Betrachtungs- 



eine Obergrenze fiir das Inkrement bei einem Fehler von err<.-^ als: a < 
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und Beleuchtungsvektoren werden im Originalvolumen ausgewertet und die 
Beleuchtungsfunktionen werden so umgeformt, dass sie im Originalvolumen 
gelten. Somit findet das Shading vollstandig im Originalvolumen statt, woraus 
sich die Bezeichnung „lnneres Shading" herleitet. 

5 

Die Schritte Shaden und Rendern sind somit gemass der vorliegenden 
Erfindung vollstandig getrennt. Das Originalvolumen v»^ird zuerst zumindest 
elementweise geshadet und dann uber Texturschnitte gerendert. Wird das 
gesamte Originalvolumen zunSchst vollstandig geshadet und dann erst 
10 gerendert, spricht man bei dieser speziellen Implementation auch von einem 
„zweistufigen Inneren Shading". 

Um das Verfahren zu implementieren, mussen die Tesselierung und die Form 
der Beleuchtungsfunktionen im Originalvolumen angeben werden. Die drei 
15 Koordinaten im Originalvolumen werden mit (a.p,y) bezeichnet und jene des 
kartesischen Projektionsraumes als (x,y,z). Die Transformation T zwischen 
den Systemen sei (x.y.z)= T(a.p,Y) - T ist l.a. kein linearer Operator, sein 
Inverses T-"" muss jedoch bis auf die Singularitaten existieren. Die 
entsprechenden linearen Operatoren werden wie folgt benannt: 

20 

X = T^(a.p.y), y = Ty(a,p.Y) und z = T^(a,p.Y). 

Im Falle von Kugelkoordinaten sind (a,p,Y) = (r,4^<p) und x = r cos((p)cos(<|>), 
y = r • cos(cp)sin(<|>) und z = r • sln((p) , fur Zylinderkoordinaten kann man 
25 (a,p.Y) = (r,<|>.z) venwenden und erhalt x = r cos{^), y =r sin(<|)) und z = z. 

Es ist aus der Literatur bekannt, dass ein Gradient in einen Belelichtungswert 
QberfQhrt werden kann, indem er als Zeiger in ein Datenfeld benutzt wird, 
welches auf den Seitenfiachen des Einheitswurfels angeordnet ist. Die Werte 
30 des Shadings werden also anhand der Vektoren der Betrachtung und der 
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Beleuchtung tabelliert und dann in einer Tabelle, Cubemap genannt. abgelegt. 
Das reduziert das Problem darauf zu zeigen, dass das Shading in 
Abhangigkeit des Gradienten im Ursprungsvolumen tabelliert werden kann. 

5 Dies ist trivialerweise gegeben - beim Fiillen der Cubemap entspricht der 
Gradientenvektor dem Veklor vom Ursprung an die zu berechnende Position 
in der Cubemap, dieser Vektor n(o) ist in Go gegeben, er kann also 

punktweise mittels T in den Projektionsraum Qberfuhrt werden: 
dn(o) 

10 

Dessen Beleuchtungsfunktion ist dann berechenbar. 



Es gilt die Definition x. :=T(xo), wobei — (x) die Richtungsableitung von T in 

an 

der Richtung von n an der Stelle x ist. Da sowohl T als auch — l a. nicht- 

dn 

15 linear sind muss die Cubemap nicht nur fur jeden Knoten der Tesselierung 
neu berechnet werden, sondern sie kann auch zwischen den Punkten nicht 
korrekt interpoliert werden. Dieser allgemeine Ansatz ist also fur krummlinige 
Koordinaten nicht gangbar. 

20 Im folgenden soli gezeigt werden, dass die Shading Werte mittels den vier 
elementaren Rechenoperationen aus den Gradienten des Ursprungsvolumens 
gewonnen werden konnen. Dies ist wegen der Ailgemeinheit der 
Transformation T nur als Naherung denkbar, die nachstehend angeben vwrd. 

25 Hierzu werden im Originalvolumen zwei Vertices Xo und x, der Tesselierung, 
zwischen denen aile Funktionen linear interpoliert werden sollen, betrachtet. 
Die kartesischen Gradienten g, := v(T(V(x,))) fur jedes i konnen per Hardware 
in Beleuchtungswerte umgerechnet werden. Es stehen jedoch nur die 
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Gradienten §01 := v(v(Xi)) zur Verfugung, die gem§ss der Erfindung nicht 
durch eine Rotation tiberfuhrt werden sollen. 

Daher mussen zunSchst die Normalen durch das Differenzial von T hergeleitet 
werden: eine Normalkomponente g in Riclitung a an der Koordinate 

=(a_p^y)wird dann im Projektionsraum durch 9 — (xo) ausgedruckt. Weil 

die Richtungsvektoren lokal definiert und differenzierbar sind, konnen sie 

linear transformiert werden, und es kann die Obergangsmatrix 

Go -^Gp(Ursprungsvolumen nach Projektionsraum) der Nonnalenvektoren im 



Punkt Xo angeben werden: N := 



da 3fi Sy 

sfT^ av, aT), 

ax dp dy 

aT, 3T, at, 

an ap 



. Dies zeigt, dass das Shading 



von nicht-karthesischen Datensatzen konventionell durchfuhrbar ist, jedoch 
erst nach der Obersetzung in den Projektionsraum. 



Daher werden die Betrachtungen im weiteren auf eine Abwandlung einer 



Basis Br •=f— ,— ,— 1 des Projektionsraumes beschrankt. 



Form der Einheitsmatrix. Die zugehorige Basis := 



In Sq hat N die 



ist 









da 










'\\ffr 






Ik 



orthonormal, wenn das Originalvolumen Go lokal orthogonal ist - was bei 
zylindrischen und spharischen Koordinatensystemen der Fall ist - und daher 
durch Drehung der Standardbasis des Projektionsvolumens Gp zu erzeugen. 



N ist in eine Diagonalmatrix mit den Elementen {q„,qp,q^} mit q» := I— 



|cT|I 




m 





Gradient in Go einem karthesischen Gradienten entspricht, welcher durch N"^ 
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transformiert wurde. Und ist die Diagonalmatrix mit den Kehrwerten dieser 
q's. Es kann also unmittelbar der Gradient in Gq venwendet werden, nachdem 
er komponentenweise mit den Kehn^^erten der q's skaliert wurde. Damit ist 
berelts eine magliche Implementation des Shadings in Go aufgezeigt. 

5 

Es sollen jedoch noch weitere Eigenschaften des Shadings ausgebeutet 
werden, so dass eine effizientere Implementierung angegeben werden kann, 
welche auch fur nicht-orthogonale Systeme gilt, welche „windschiefe" 
Basisvektoren besitzen. Solche Systeme konnen beispielweise von 
10 sogenannten Array Scannem erzeugt werden. Da im Stand der Technik nie 
die Gradienten selbst verwendet werden, sondern alle Beleuchtungswerte aus 
den Skalarprodukten (g,!), (g,v) und (v,T) und einer optionalen Skalierung mit 
||g|| erzeugt werden, ergeben sich hier Moglichkeiten zur Eihsparung. 

15 Hierzu mussen die Gleichungen lediglich im Originalvolumen in eine 
geeignete Form gebracht werden. Die einfache Forderung lautet dann 
(g. j)= (goi. b) (^^- 3) etc., also die Invarianz der Skalarprodukte. Diese 
Forderung ist trivialerwelse erfullbar, indem dasselbe Euklldische 
Skalarprodukt verwendet wird, welches basisunabhangig ist. Dies ist aber kein 

20 gangbarer Weg, well ein solches Skalarprodukt aufwendig zu berechnen ist. 
Das Kriterium fur den Ldsungsansatz ist einfach: g soil nicht transformiert 
werden, well fur jeden Pixel ein g anfallt, g und v werden jedoch nur in jedem 

Vertex ausgewertet. Es mussen also i^ij)p = (gohS(t))o (ep- 4) zur Verfugung 
gestellt werden, mit einem basisabhangigen Standardskalarprodukt fur das 
25 Originalvolumen und den Projektionsraum. 

Im weiteren wird X •= s(k) explizit hergeleitet. 
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Im betrachteten Punkt x im Ursprungsvolumen Gq seien die Koordinaten der 
Vektoren durch die Basis e:= ^„,ep,ej gegeben. also g:= g.e„ +gpep +gye^ 
und T := l„e„ + ipOp +l^e^ . Betrachtet man die Skalarprodukte der Basis 
Cg := (iilij), dann ist Ch = 1 . aber i.a c, 0 fur die nicht-ortliogonalen Basen. 
5 Ausgeschrieben lauten die Skalarprodukte: 

(ftT)= gja +gp'p +9t'y ^-gaipc^p +gaiyC^ +gpi«c«p +gpiyCp, +g,iaCo^ +gY'pCpr (eq- 5) 

wobei bei orthonormaler Basis e nur die ersten drei Terme existierten. 

10 

Gruppierl man den Ausdruck in (eq. 6) nach den Komponenten von g: 

(g.T)= ga(la +lpCap +lrc„,)+gp(lp +i„c^ + l,Cp,)+ g,(l, -f-l^c^ +lpCp,) (eq. 6), 

15 erkennt man (eq 6) als Standardskalarprodukt mit einem neuen Vektor 1' mit 
den Komponenten 

S(l) = (l„ + IpC^ + l,c^,lp + Uc^ + l^cp,.!^ + l„c^ + IpCp,) (eq. 7). 

20 Mit S(l) steht also eine universale Transformationsvorschrifl fur die Vektoren I 
und V zur Verfugung. Diesen neuen I' und v' darf man keine wirkliche 
geometrische Bedeutung zuschreiben, ausser dass sie die Skalarprodukte im 
Originalvolumen Gq invariant erhalten. 

25 Damit steht eine sehr einfache Vorschrift von minimalem Rechenaufwand zur 
Verfugung, welche das Berechnen der Skalarprodukte im Originalvolumen Go 
auf dieselbe Art wie im karthesischen Projektionsraum eriaubt, ohne dass T 
hierzu ausgewertet werden musste. Die „c"s sind ihrerseits Skalarprodukte der 
lokalen Basen, werden fiir die meisten Systeme - auch fur nicht orthogonale - 
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jedoch globale Konstanten sein, und daher keinerlei Rechenaufwand 
konsumieren. 



I und V warden in und und dazwischen vom Grafikprozessor linear 
5 interpoliert, die „g"s hingegen werden in jedem Pixel berechnet und nicht nur 
In den Vertices. Der durch die Linearislerung eingefuhrte Fehier existiert beim 
kartesischen Rendem ebenfalls, da der Sehvektor uber das Sehfeld 
schwenkt, das Schwenken aber als lineare Operation angenommen wird. 

10 Dennoch ist der Fehier der Interpolation, der allgemeinen Abbildung T wegen, 
viel komplexer als im karthesischen Fall: der Beleuchtungsvektor I und der 
Betrachtervektor v sind beides Richtungsvektoren, welche durch Normierung 

in einem Punkt x erzeugt werden v := .^"^^ (eq. 8) mit x^ der 

Ir ~ II 

Betrachterposition. Um die einfachst mogliche Form der Richtungsvektoren zu 
15 erhalten, verwenden wir (eq. 7). 

Somit herrschen In Go dieselben Bedingungen wie in karthesischen 
Systemen. Da (eq 3) in Gq linear ist, kommt die Nichtlinearitat einzig vom ' 
Nenner von (eq. 8). Fur eine „hinreichend feine" Tesselierung in dem Sinne, 
20 dass sich der Nenner von (eq. 8) monoton entwickelt, kann der Fehier der 
Linearsierung durch die Kehrwertfunktion abgeschatzt werden: 



e(t): 



°4'i^'^-'^i]-[ ..d„.(i-.).d. ] f^-^'- '^t°-^i '^"^ 

Interpolationsparameter und die „d"s stehen fur die Nenner in den Vertices, 
und Q ist eine Skalierung, welche daher kommt, dass die durch (eq. 7) 
25 umgeformten Vektoren keine Einheitsvektoren mehr sind. Durch die 

Beschrankungen der „c"s in a3 gilt Q g ]d..4[. Da e(t) = O.t e {0,1}, kann das 
Fehlermaximum an den Extreme von e ausgewertet werden: 
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^^O^J 1+ ^0 = 0 (eq. 10a), was fur den nicht trivalen 

dt do d, (t.do+(l-t).dO' 

Fall, dass do ^ d^ (andernfalls ist der Fehler identisch 0) ausgewertet werden 
zu: 

5 (t(do-d,)+dJ=-AlA = ^zd2.dod,=dod, unddamit 

d^"d; 



do -d, 



(eq. 10b). 



Die Losung Ist eindeutig, da t e [0..1] und daher ein Pluszeichen vor der 
10 Wurzel stehen muss. 

Dies ergibt ein handliches Kriterium, die Feinheit der Tesselierung 
festzulegen. 

15 Durch Ruckelnsetzen von (eq. 10b) in (eq. 1 0a) erhalt man: 



e(t):=Q 



Icio dj d, 



do-di 



= Q 



rVdodi-d,] 


1 1 


= Q 


J 1_ 


t dod, J 


^ d, ^dodi 




di do 



20 Das ergibt das uberraschende Resultat, dass ahnlich wie fur (eq 1 ), Fehler<1 
erzwingbar ist. Wenn do vorgegeben ist, so bestimmt (eq 11), in welchem 
Berelch d, gewahit werden kann. 
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Es ist sogar so, dass mit der Tesselieaing verhaltnismassig sorglos 
umgegangen werden kann, solange man sich im Originalvolumen Go bewegt, 
well man dann keiner Singularitat nahekommt und die „d"s verhaltnismdssig 
gross blerben. Fur „di''s grosser als 1 , ist die Methode durch den 
5 Geometrischen Fehler von (eq 1) begrenzt und nicht durch die Fehler der 
Shading Naherung (eq 9). 

Bisher wurde das folgende gezeigt: die fOr das allgemeine Shading und 
Rendering notwendigen Skalierungen und Transformationen sind mit 

10 elementaren Rechenoperationen durchfuhrbar. In den Vertices erhalt man die 
korrekten Beleuchtungswerte, die interpolierten Werte enthalten hochstens 
Fehler zweiter Ordnung wie in (eq 9) beschrieben. Die Auswertung in den 
Vertices ist ebenfalls elementar, und die Vertices sind vom Standpunkt des 
Rechenaufwandes nicht sehr dicht (und mQssen auch nicht sehr dicht sein); 

15 ihre Auswertung ist folglich nicht der Flaschenhals der Berechnung. 

Im folgenden wird nun die Erfindung anhand der schematischen Zeichnung 
ndher erldutert Es zeigen: 

t 

J 

20 Fig. 1 Ein aus dem Stand der Technik bekanntes Verfahren zum 

Visualisierien eines Originalvolumes; 

Fig. 2 Erzeugung einer ebenen Textur gemass dem Stand der Technik; 

Fig. 3 Linearisierung einer kugelschalenformigen Textur; 

Fig. 4 ein Ausfuhrungsbeispiel des erfindungsgemassen Verfahrens; 



25 



Fig. 5 



Beispiele fur Texturen in Zylinderkoordinaten; 
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Fig. 6 Ausftihrungsbeispiel fur eine Linearisierung in krummlinigen 
Koordinaten. 

Die Figuren 1 bis 3 zeigen den Stand der Technik und wurden bereits 
Eingangs eingehend diskutiert. 

Fig. 4 zeigt schematisch die wichtigsten Scliritte eines Ausfulirungsbeispiels 
des erfindungsgemassen Verfalirens. Anhand von Fig. 4 soil beispielliaft auf 
den fur die Praxis wichtigen Fall Bezug genommen werden, dass ein 
Datensatz D auf Messwerten D(a, p. y) beruht, die aus einer 
volumenaufgeiesten Abtastung eines Korpers resultieren, und dass aus 
diesem Datensatz D eine dreidimensionale Darstellung erzeugt wird. Mit 
"dreidimensionaler Darstellung " ist dabei gemeint, dass die Darstellung 
tatsachlich dreidimensional ist, oder dass beispielsweise durch eine geeignete 
Projektionsvorrichtung eine stereoskopische Projektion erzeugt wird, oder 
dass die Projektion eben, beispielsweise auf einem Computermonitor erfolgt, 
aber ein dreidimensionaler Eindruck vermittelt wird, z.B. mittels Methoden der 
rSumlichen oder perspektivischen Darstellung, insbesondere unter 
Verwendung eines entspreclienden Beleuchtungsmodells. Im speziellen 
konnen solche Darstellungen semi-transparent sein, so dass sle einen 
Einblick in den abgetasteten Korper gestatten. 

Bei dem in Fig. 4 dargestellten AusfOhrungsbeispiel wird ein Originalvolumen 
Go, beispielsweise das Herz Go eines menschlichen KGrpers, mit einer 
Messeinrichtung 1 , im vorliegenden Beispiel mit einer Ultraschall- 
Messeinrichtung 1 , abgetastet. Eine solche Ultraschall-Messeinrichtung 1 
umfasst beispielsweise mehrere Ultraschallwandler 1 1 , von denen mehrere 
bezuglich einer Achse y benachbart zueinander angeordnet sind. Wahrend 
des Betriebs wird die Ultraschall-Messeinrichtung 1 um die Achse y rotiert, wie 
dies der Doppelpfeil an der Achse y der Ultraschall-Messeinrichtung 1 in Fig. 4 
andeutet. Dabei werden die einzelnen Ultraschallwandler 1 1 im wesentlichen 
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parallel betrieben, so dass das Originalvolumen Go, also beispielsweise ein 
Herz, gleichzeitig in mehreren parallelen Schichten, die jeweils im 
wesentlichen senkrecht zur y - Achse liegen, uber einen Winkelbereich oder 
einen Sektor (Po X yo) mit Ultraschall abgetastet wird. Die Information 
5 bezuglich der diitten Dimension in Richtung a wird beispielsweise aus der 
Laufzeit der Ultraschailechos gewonnen. 

Die so gewonnenen volumenaufgelosten Daten D(a, p, y) liegen, 
gegebenenfalls nach elner signaltechnischen Vorverarbeitung, somit als ein 
raumlich aufgeloster Datensatz D vor, der die Informationen der 

10 abzubildenden Struktur enthalt. Dabei ist der Datensatz D durch sognannte 
Texturen Tp aufgebaut, die Schnitte durch das Orginalvolumen Go 
reprasentieren. Durch die Ultraschall-Messeinrichtung 1 selbst bzw. durch 
deren Arbeitsweise ist das Messkoordinatensystem Km, in welchem die Daten 
D(a, p, y) vorliegen, vorgegeben. Bei dem hier beschriebenen Beispiel handelt 

15 es sich um Zylinderkoordinaten. Die Texturen Tp entsprechen dann, wie in* 
Fig. 5 gezeigt, drel verschiedenen Flachentypen, also Textur-Typen Ta, Tp 
und Ty, auf denen jeweils eine der Zylinderkoordinaten a, p, y einen 
konstanten Wert hat Dabei ist die Anwendung des erfindungsgemassen 
Verfahrens selbstverstandlich keineswegs auf Zylinderkoordinaten 

20 beschrankt, sondem auch auf andere krummlinige Koordinaten und 

selbstverstandlich auch auf karthesische Koordinaten anwendbar und das 
zuvor gesagte gilt fur andere als Zylinderkoordinaten vollig analog. Der durch 
die Ultraschall-Messeinrichtung 1 generierte Datensatz D wird in eine 
Datenverarbeitungsanlage 3 geladen und dort verarbeitet Dabei wird, 

25 eventuell nach einer mathematischen Aufbereitungsphase, zunachst das 
Beleuchtungsmodetl, d.h. die diesem zugrunde liegenden 
Beleuchtungsfunktionen im Messkoordinatensystem Km ausgewertet, d.h. es 
wird das oben definierte "Innere Shading" durchgefuhrt. 
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Danach werden die Daten dem ebenfalls oben definierten "Rendern" 
unterzogen und der Grafikhardware 4 zugefuhrt, mit deren Hilfe eine 
dreidimensionale Darstellung des zu projizierenden Bildes 5, beispielsweise 
auf einem Beobachtungsmonitor 2 erfolgt. 

5 Fig. 5 zeigt schematisch die drei Textur-Typen, fCir den Fall , dass das 
Messkoordinatensystem Km Zylindersymmetrie besitzt. Besitzt das 
Messkoordinatensystem Km eine andere Symmetrie, beispielsweise 

Kugelsymmetrie, haben die Texturen Ta, Tp und Jy bzw. die sie 

reprasentierenden Flachen selbstverstandlich entsprechend der Symmetrie 

10 des Koordinatensystems, eine andere Geometrie. 

Erfindungsgemass wird nach dem Durchfuhren des Shadings, d.h. nach dem 
Auswerten des Beleucintungsmodells BM im Messkoordinatensystem Km, der 
Datensatz D, bzw. werden die Texturen Tp im Messkoordinatensystem Km 
linearisiert, d.h. dem Vorgang des Renderns unterzogen. Dies soil beispielhaft 
15 in Fig. 6 in einer schematischen Darstellung an einer Gruppe von Ta Texturen 
eines Datensatzes in einem zylindersymmetrischen Messkoordinatensystem 
Km noch kurz eriautert werden. 

Wie bereits eingehend eriautert, wird der Sehvektor S (bzw. der 
Beleuchtungsvektor) weder auf ebenen Texturen Tp noch auf krummlinig 
20 begrenzten Texturen Tp einen konstanten Betrag haben; vielmehr ist der 

Betrag dieser Vektoren, wie im linken Bild von Fig. 6 schematisch dargestellt, 
auf Kugelschalen K konstant, Somit andert sich der Sehvektor S im 
allgemeinen uber einen Schnitt, d.h. uber eine Textur Tp. 

Well somit betragsmassige Anderungen des Sehvektors S in der Regel nicht 
25 linear von den Texturkoordinaten der Texturen Tp, abhangen, konnen diese 
auch nicht auf der Graphikhardware 4 unmittelbar kodiert werden. Daher 
werden die krummlinig begrenzten Texturen Tp, wie im mittleren Bild von 
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Fig. 6 dargestellt, linearisiert, so dass geradlinig begrenzte Teilflachen 
entstehen. Die Eckpunkte dieser Teilflachen konnen nun ohne eine 
Koordinatentransfomnation durchzufuhren auf der Graphikhardware 4, wis im 
rechten Bild von Fig. 6 schematisch dargestellt, kodiert und beispielsweise auf 
5 einem Beobachtungsmonitor 2 als Bild 5 dargestellt werden. 

Durch das erfindungsgemasse Verfahren ist es somit moglich, auf besonders 
elegante und effiziente Weise mehrdimensionale Datensatze, die in beiiebigen 
krummlinigen Koordinaten, z.B. in Zylinder- Oder Kugelsymmetrie vorliegen, 
ohne grossen Rechenaufwand und enorm schnell auf gewohnlicher 

10 Graphikhardware unter Berucksichtigung eines entsprechenden 

Beleuchtungsmodells unmittelbar zu kodieren. Das wird dadurch erreicht, 
dass der Vorgang des sogenannten Shadings, das heisst der Auswertung der 
Beleuchtungsfunktionen und des Rendems, d.h. der Vorgang der 
Linearisierungsoperationen, vollstSndig getrennt sind: Die 

15 Beleuchtungsfunktionen werden im Originalvolumen, also im 

Messkoordinatensystem ausgewertet und dann erst linearisiert. Da die 
Auswertung der Datensatze m it elementaren Rechenoperationen moglich ist 
und das Shading ausschliesslich im Originalvolumen und nicht erst, wie aus 
dem Stand der Technik bekannt, im Projektionsraum stattfindet, ist das 

20 erfindungsgemasse Verfahren enorm schnell, so dass in Echtzeit die 

Visualisierung auch von sehr komplexen Datensatzen m6glich wird. Selbst die 
Visualisierung bewegter Vorgange in stereoskopischer Darstellung wird mit 
hoher Auflosung und mit Bildwiederholraten, die typischen Videofrequenzen 
entsprechen, moglich. Dadurch, dass die Beleuchtungsfunktionen im 

25 Originalvolumen und nicht im Projektionsraum ausgewertet werden, ist es mit 
dem erfindungsgemassen Verfahren aufgrund seiner beachtllchen 
Geschwindigkeit bei der Durchfuhrung der Visualisierungsoperationen 
moglich, auch sehr komplexe Beleuchtungsmodelle zu verwenden, so dass 
hochaufgeloste und realistische Darstellungen von bisher nicht gekannter 

30 Qualitat erreichbar sind. 
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Patentanspruche P.7321/lr/uw 

1 . Verfahren zum Visualisieren eines raumlich aufgelosten Datensatzes (D) 
unter Verwendung eines Beleuchtungsmodells (BM), wobei ein Datum 
(D(a,P,Y)) des Datensatzes (D) jeweils einem Volumenelement (V) 

5 zugeordnet ist, dessen Lage durch Koordinaten (a, p, y) in einem 
Messkoordinatensystem (Km) beschrieben wird, wobei die Daten 
(D(a,p,Y)) als mindestens eine Textur (Tai, Tpj, Tyk) in eine Grafikhardware 
geladen werden, um in einem Projektionsraum eine bildliclie Darstellung 
(5) zu erzeugen, dadurch gekennzeichnet, dass das Beleuchtungsmodell 

10 (BM) im Messkoordinatensystem (Km) ausgewertet wird. 

2. Verfahren nach Anspruch 1 , bei welchem die Daten (D(a,p,y)) des 
Datensatzes (D) ohne Transformation aus dem Messkoordinatensystem 
(Km) in ein anderes Koordinatensystem, insbesondere ohne 
Transformation in ein karthesisches und / oder isotropes 

15 Koordinatensystem, verarbeitet werden. 

3. Verfahren nach einem der vorangehenden Anspruche, bei welchem das 
Messkoordinatensystem (Km) ein nicht-karthesisches 
Messkoordinatensystem (Km) ist. 

4. Verfahren nach einem der vorangehenden Anspruche, bei welchem das 
20 Messkoordinatensystem (Km) ein Zylinder- oder Kugelkoordlnatensystem 

(Km) ist. 



5. 



Verfahren nach einem der vorangehenden Anspruche, bei welchem 
zwischen den Daten (D(a,p,y)) des Datensatzes (D) im 
Messkoordinatensystem (Km) linear interpoliert wird. 
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6. Verfahren nach einem der vorangehenden Anspruche, bei welchem das 
Beleuchtungsmodell im Datensatz (D) nahe einer Singularitat ausgewertet 
wird. 

7. Verfahren nach einem der vorangehenden Anspruche, bei welchem die 
5 Daten (D(a,p.Y)) des Datensatzes (D) eine volumenaufgeloste Abtastung 

eines Korpers (Go) reprasentieren, und bei welchem die bildliche 
Darstellung (5) eine dreidimensionale, insbesondere eine semi- 
. transparente Darstellung (5) des Korpers (Go) ist. ; v. 

8. Verfahren nach einem der vorangehenden Anspruche, bei welchem die 
10 bildliche Darstellung (5) als stereoskopische Projektion erzeugt wird. 

9. Verfahren nach einem der vorangehenden Anspruche, bei welchem die 
Daten D((a,p,y)) des Datensatzes (D) mittels einer Ultraschall- 
Messeinrichtung (1) generiert werden. 

10. Verwendung eines Verfahrens gemass einem der vorangehenden 
15 Anspruche, insbesondere fur medizinische Zwecke, zum schnellen 

Generieren dreidimensionaler Darstellungen (5) eines Korpers (Go), 
insbesondere eines menschlichen Korpers oder Teilen davon, anhand von 
messtechnisch gewonnenen Daten (D(a,p,Y)). 



20 
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Zusammenfassunq P.7321/lr/uw 

Erfindungsgemass wird ein Verfahren zum Visualisieren eines raumlich 
aufgelosten Datensatzes (D) unter Verwendung eines Beleuchtungsmodells 
(BM) vorgeschlagen, wobei ein Datum (D(a.p,y)) des Datensatzes (D) jeweils 

5 einem Volumenelement (V) zugeordnet ist. dessen Lage durch Koordinaten 
(a, p, y) in einem Messkoordinatensystem (Km) beschrieben wird. Die Daten 
(D(a,p,y)) warden als mindestens eine Textur (Toi, Tft, Tyk) in eine 
Grafikhardware (4) geladen, um in einem Projektionsraum eine bildliche 
Darstellung (5) zu erzeugen. Dabei wird das Beleuchtungsmodell (BM) im 

10 Messkoordinatensystem (Km) ausgewertet. 



(Fig: 4) 
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